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We study the charge transport properties of triangular quantum dot molecule (TQDM) connected to metallic elec¬ 
trodes, taking into account all correlation functions and relevant charging states. The quantum interference (QI) 
effect of TQDM resulting from electron coherent tunneling between quantum dots is revealed and well interpreted by 
the long distance coherent tunneling mechanism. The spectra of electrical conductance of TQDM with charge filling 
from one to six electrons clearly depict the many-body and topological effects. The calculated charge stability dia¬ 
gram for conductance and total occupation numbers match well with the recent experimental measurements. We also 
demonstrate that the destructive QI effect on the tunneling current of TQDM is robust with respect to temperature 
variation, making the single electron QI transistor feasible at higher temperatures. 


1 Introduction 

Molecule transistors (MTs) provide a brightened scenario 
of nanoelectronics with low power consumption. 1 To 
date, the implementation of MTs remains challenging, 
and a good theoretical understanding of their charac¬ 
teristics is essential for advancing the technology. The 
current-voltage (I-V) curves of MTs are typically pre¬ 
dicted by calculations based on the density functional 
theory (DFT). 3 However, the DFT approach cannot fully 
capture the correlation effect in the transport behavior 
of MTs in the Coulomb-blockade regime. A theoretical 
framework to treat adequately the many-body problem 
of a molecular junction remains elusive due to the com¬ 
plicated quantum nature of such devices. Experimental 
studies of a artificial molecule with simplified structures 
are important not only for the advances of novel nanoelec¬ 
tronics, but also for providing a testing ground of many 
body theory. For example, the coherent tunneling be¬ 
tween serially coupled double quantum dots (DQDs) was 
studied and demonstrated for application as a spin filter 
in the Pauli spin blockade regime. 4 Recent experimen- 
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tal studies have been extended to serially coupled triple 
quantum dots (SCTQDs) for studying the effect of long 
distance coherent tunneling (LDCT) in electron trans¬ 
port] 5 7 Triangular quantum dot molecule (TQDM) pro¬ 
vides the simplest topological structure with quantum 
interference (QI) phenomena.^®^ The QI effect in the 
coherent tunneling pro cess o f TQDM junctions has been 
studied experimentally.^^ 

It was suggested that the tunneling currents through 
benzene molecules can also show a destructive QI behav¬ 
ior. 39 The tunneling current through a single benzene 
molecule was theoretically studied by DFT. 9 However, 
the influence of the strong correlation on the QI effect 
remains unclear due to the limitation of DFT. Many the¬ 
oretical works have pointed out that electron Coulomb 
interactions have strong influence not only on the elec¬ 
tronic structures of TQDM, 13 14 but also on the prob¬ 
ability weights of electron transport paths. 15 18 When 
both the intradot and interdot Coulomb interactions in 
a TQDM are included, the transport behavior involving 
multiple electrons becomes quite complicate. The setup 
for the TQDM junction of interest is depicted in Fig. 1. 
Here we present a full many-body solution to the tunnel¬ 
ing current of TQDM, which can well illustrate the Pauli 
spin blockade effect of DQDs, 4 LDCT of SCTQDs 5-7 
and QI of TQDMs 11 for both equilibrium and nonequilib¬ 
rium cases. Thus, our theoretical work can provide useful 
guidelines for the design of future molecular electronics 
and the realization of large scale quantum registers built 
by multiple QDs.^ 

We adopt the equation of motion method (EOM), 
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Hubbard Hamiltonian for multiple QDs. 
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Fig. 1 Illustration of the TQDM junction system of 
interest. 

which is a powerful tool for studying electron transport, 
taking into account electron Coulomb interactions. 15 18 
This method has been applied to reveal the transport 
behaviors of a s ingle QD with multiple energy levels^ 
and DQDs.E^ 18 For a TQDM with one level per QD, 
there are 64 configurations for electrons to transport be¬ 
tween electrodes] 20 Previous theoretical works have ig¬ 
nored the high-order Green functions resulting from elec¬ 
tron Coulomb interactions to simplify the calculation. 20 
To have a full solution becomes crucial for depicting the 
charge transport involving a few electrons. We solve the 
EOM of Green functions up to six electrons, taking into 
account all correlations caused by electron Coulomb in¬ 
teractions and electron hopping between TQDM. This 
involves solving 4752 Green’s functions and 923 correla¬ 
tion functions self-consistently. 

2 Model 

We consider an artificial molecule made of nanoscale 
QDs, in which the energy level separations are much 
larger than the on-site Coulomb interactions and ther¬ 
mal energies. Thus, only one energy level for each quan¬ 
tum dot is included. The extended Hubbard-Anderson 
model is employed to simulate the TQDM junction with 
Hamiltonian given by H = Hq + Ht + Hqd s , where 
Ho = J2k,a,a e k c i,a,a c k,cr,a is the Hamiltonian for free 
electrons in the electrodes. Cj[ (ck,o-,a) creates (de¬ 
stroys) an electron of momentum k and spin a with en¬ 
ergy e k in the a electrode. H T = , a de,<r + 

KaA, a c k ,cr,a)- Vk,a,£ describes the coupling between 
the a electrode and the l-th QD. d\ a (dq a ) creates (de¬ 
stroys) an electron in the Cth dot. Hqds is the extended 


H QDs ^ ^ r T ^ ^ (l) 

£,<j i 

£<j,cr,cr' 

where Eg is the spin-independent QD energy level, = 
d\ a d J £, (7 , Ui and Uij (£ < j ) denote the intradot and 
interdot Coulomb interactions, respectively and tij de¬ 
scribes the electron interdot coupling. The interdot 
Coulomb interactions as well as intradot Coulomb inter¬ 
actions are important for nanoscale semiconductor QDs 
and molecules. Therefore, Uij cannot be ignored. 

Using the Keldysh-Green’s function technique^ 5 2 -^, the 
electrical current from reservoir a to the TQDM junction 
is calculated according to the Meir-Wingreen formula 

J a = j J de'£r?[Gf tr {e) + f a (e)(G r ja (e) - G^(e))}, 

jo 

( 2 ) 

where T“(e) = \Vk,a,j\ 2 d( e — e k ) is the tunneling rate 
between the <a-th reservoir and the j-th QD. Through¬ 
out the paper, for two-terminal devices we assume that 
the left (right) lead is only coupled to the left (right) 
QD with tunneling rate T l (T#), while there is no cou¬ 
pling between the center QD and the two leads. For 
three-terminal devices, the coupling between the cen¬ 
ter QD and a third gate is described by the tunnel- 
ing rate T c . / Q (e) = l/{exp[(e - l^ a )/k B T] + 1} de- 
notes the Fermi distribution function for the <a-th elec¬ 
trode, where (i a is the chemical potential and T is 
the temperature of the system. e, ft, and ft# de¬ 
note the electron charge, the Planck’s constant, and 
the Boltzmann constant, respectively. G^(e), GJ a (e), 
and Gj a (e) are the frequency domain representations of 
the one-particle lessor, retarded, and advanced Green’s 
functions Gf a (t,t') = i{d\ (J {t l )d^ (J {t)), G r ja (t,t') = 

-iO(t - and G a jrj {t,t’) = iO(t'~ 

t)({dj )(7 (t),dj a (t')}), respectively. These one-particle 
Green’s functions are related recursively to other Green’s 
functions and correlators via the many-body equation 
of motion,which we solve via an iterative numer¬ 
ical procedure to obtain all n-particle Green’s func¬ 
tions (n = 1, • • • ,6) and correlators for the TQDM. 
(See supplemental materials.) Our procedure is valid 
in the Cou lomb blockade regime, but not the Kondo 
regime. 22 23 Throughout this paper, we assume the on¬ 
site Coulomb interaction Ug = Uo = lOOTo for all three 
QDs and the same tunneling rates at all leads, T™ = T 
with j labeling the QD directly connected to lead a. 
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Fig. 2 Electrical conductance (G e ) of TQDM as a function 
of central QD energy for different t^R strengths in the Pauli 
spin blockade configuration with El — E F and 
E r = Ef - Ur. (a) the = tcR = 3r 0 , and (b) tgj replaced 
by —tgj. Other physical parameters are k F T = ir 0 , 
r — 0.3r o , Ulc — Ucr — 30r o , and Ulr — 0. 


3 Results and discussion 

Although the QI of TQDM was theoretically investi¬ 
gated previously 24 ^, the effect of intradot and inter¬ 
dot Coulomb interactions was not considered. Here, we 
utilize the LDCT effect to tune the effective hopping 
strength between outer QDs to achieve the destructive 
and constructive QIs in the presence of electron Coulomb 
interaction. We consider a TQDM junction in the Pauli 
spin blockade (PSB) configuration 4 181 with E l = E f , 
Er = E f — Uo, tLc = tcR = 3T 0 , Tl = Tr = 0.3T 0 , 
Ulc = Ucr = 30T 0 , and U LR = 0 . Fig. 2 (a) shows 
the electrical conductance (G e ) as a function of central 
QD energy level (A q = Ec — E f ) for t lr varying from 
0 to 0 . 2 T 0 at kpT = 1 T 0 (in weak interdot coupling 
regime). For A c less than — 15Fo, G e is not sensitive 
to the variation of indicating that the transport is 
mainly through the upper path involving the center QD 
as shown in the inset of Fig. 2(a). In the case of t R R * 0, 
G e can be well explained by the LDCT effect when Ec is 
far away from E f The central QD provides an interme¬ 
diated state for electrons in the outer QDs. Through the 
upper path, TQDM behaves like a double QD with an ef¬ 
fective hopping t e ff = -tLctcR/(U C R + Ac), which can 
also be understood by the second order perturbation the¬ 
orem] Once tLR 7 ^ 0 (the lower path turns on), electron 
transport through the two paths with t e ff and tLR lead 
to a destructive QI. Note that G e for the case of tLR 7 ^ 0 
is reduced compared to the case of tLR = 0. In particular, 


G e is vanishingly small at the value of Ac where \t e ff\ 
crosses tLR as indicated by dashed lines in Fig. 2(a). For 
illustration, the curve for \t e ff\ = \t L ctcR\/{UcR + Ac) 
is also shown in Fig. 2 (a). (See short-dashed curve) 
The vanishing G e occurs at lower Ac with increasing 
| tLR\ (Compare dash-dotted with dashed curves). Due 
to topological effect, the electron-hole symmetry does not 
hold for the energy spectrum of TQDM . 13 When TQDM 
has identical QD energy levels {Eg = Eq) and homoge¬ 
nous electron hopping strengths tgj = £, we have one 
level 6 = Eq T 2£ and one doubly degenerate level with 
e = Eq — t for the case of Ug = Ugj = 0. Unlike the cases 
of DQDs and SCTQDs, the lowest energy level depends 
on the sign of tgj . This is a manifesting result of electron- 
hole asymmetrical behavior of TQDM. Therefore, it will 
be interesting to examine the sign effect of tgj on the QI 
behavior. Physically, the sign of tgj depends on the sym¬ 
metry properties of orbitals, which can change in different 
configurations. We can replace tgj by —tgj to examine 
the QI effect with respect to the electron-hole symmetry. 
The results are shown in Fig. 2(b). We find that G e is 
enhanced with increasing tLR, which is attributed to the 
constructive QI effect, in contrast to the destructive QI 
effect shown in Fig. 2(a). 

To gain deeper understanding of the destructive and 
constructive QI shown in Fig. 2 , we compare our 
full calc ulation with the weak interdot-coupling the- 
ory HSEH w hich allows simple closed-form expression for 
the electrical conductance of TQDM. We obtain G e = 

2 e 2 /hfdeT(e)[df{e)/dE F \ « (2 e 2 /h)T{E F ) at low- 

temperature limit, where the transmission coefficient 
T(e) with 64 configurations is approximately given by 

Tpsb(U) (3) 

__ 4T lFrPpsbEqi _ 

|MlM2M3 - t 2 CR H 1 - t 2 LC li 3 - t 2 LR H 2 - 2t LR t LC t C R\ 2 ’ 

where Fqi = l^i^LctcR,/ 1^2 + t lr ) 2 is a factor related 
to QI. hi = e - E l + iT L , H 2 = e - E c - Urc and 
H 3 = e — Er — Ur + ITr. Ppsb denotes the probability 
weight in the PSB configuration.^ From Eq. (3), we have 

^ _ PpsB^E 2 {t e ff + tLR ) 2 _ 

h (r 2 + 2 tefftLR + t 2 LR) 2 + r 2 (thl + U 2) 2 ’ 

(4) 

where t eff = -t LC tcR/(U C R+Ac),thi = -t 2 L c/(UcR + 
A c ), and t h2 = -t 2 CR /{U C R + A c ) with Tr = T# = T. 
For t lc = tcR = t c = 0, Eq. (4) reduces to the con¬ 
ductance of DQD , 17 while for tLR = 0, it reduces to the 
G e of SCTQDP 6 At Ac = 15Fo and 60Fo, which sat¬ 
isfy the condition of t e ff + t R R = 0 for tLR = 0 . 2 TO 
and tLR = OTTo, respectively, and we see G e vanishes 
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there. This well illustrates the destructive QI seen in 
Fig. 2(a). Once we make the substitution tgj —tgj 
in Eq. (4), we can reveal the constructive QI in G e as 
shown in Fig. 2(b). Note that in the weak coupling regime 
(t e ff/T < 1), the probability weight Ppsb of Eq. (4) cal¬ 
culated according to the procedures in Ref. [18j where the 
interdot two-particle correlation functions are factorized 
as the product of single occupation numbers, is consis¬ 
tent with the full calculation, but not for t e ff/T 1. 
Away from the weak coupling regime, the inter dot elec¬ 
tron correlations become important. To explicitly reveal 
the importance of electron correlation effects, we plot the 
curve of tpR = 0.2Fo (with triangle marks) calculated by 
the procedure of Ref. 18 (including 64 configurations) 
in Fig. 2(a). Comparison between the full solution and 
the approximation considered in Ref. 18, we find that 
the electron correlation effects become very crucial when 
t e ff/T 1. Once electron transport involves more elec¬ 
trons, the high-order (beyond two-particle) Green func¬ 
tions and correlation functions should be included (see 
the results of Fig. 3). The difference between the con¬ 
ventional mean-field theory of Ref. 3 with the full solution 
is even larger. The comparison between mean-field the¬ 
ory and the procedure of Ref. 18 has been discussed in 
the appendix of Ref. 18. 

According to Eq. (4), constructive and destructive QI 
effects depend on the sign of tgj. Therefore, if the wave- 
function of the center dot has opposite parity (say, an 
x-like state) with respect to the wavefunctions in two 
outer dots (assumed to be 5-like), then tpc will have op¬ 
posite sign compared with tcR and and the sign of 
t e ff will be flipped. Consequently, the destructive QI 
shown in Fig. 2(a) will become constructive QI. Thus, 
for a center QD with an 5-like ground state and an ex¬ 
cited p-like state, it is possible to see the change of QI 
between destructive and constructive by tuning the gate 
voltage, which sweeps through different resonance ener¬ 
gies of the center QD in addition to the change of sign of 
t e f f when the Fermi level goes from below the resonance 
level to above the resonance level. From the results of 
Fig. 2, QI effect can be electrically controlled by the en¬ 
ergy level E c . This advantage of TQDM may be useful 
for improving the spin filtering of DQDs®. 

In addition to QI effects, spin frustration and topo¬ 
logical effects (due to electron-hole asymmetry) on the 
measured quantities (electrical conductance or current) 
are also interesting issues.^ 1 Fig. 3(a) shows G e as a 
function of gate voltage, for Eg = Ep + lOTo — A^; 
£ = L, C, R at two different temperatures, kpT = 0 and 
lTo- Here, we consider the homogenous configuration 
with U LR = U LC = U C r = 30r 0 and tgj = t c = 3T 0 . At 
low temperature, there are six main peaks in the G e spec- 



Fig. 3 (a) Electrical conductance of TQDMs as a function 
of gate voltage A g with El «= Er = Ec — Er + lCTo — A g 
and T = 1.6ro.(b) Correlation functions for ksT m. CTo. 
Other physical parameters are tgj = 3Fo and 
Ulr — Ulc — Ucr — 30r o . 

trum, labeled by e n , n = 1, • • • ,6 and some secondary 
peaks. At higher temperature, the six main peaks are 
suppressed and broadened as shown by the dashed curve, 
and the secondary peaks are washed out. We also plot the 
total occupation number, N = ^2 a {NL, a + Nr, a + Nc, a ) 
as the black solid curve, which shows a stair-case behav¬ 
ior with plateaus at N = 1, • • • ,6, corresponding to the 
filling of TQDM with 1 to 6 electrons. It can be seen that 
the six main peaks occur at A g where N is increased by 1. 
Thus, the peak positions e n correspond to the chemical 
potential of electrons in TQDM, i.e. the energy needed 
to add an electron to the system. The main peak posi¬ 
tions can be approximately obtained by the calculation 
of chemical potential of TQDM without considering the 
coupling with leads as done in Ref. [HD For example, 
ei = El — 2|£ c |, 62 = Er T Ulc — 8t 2 /(Uo — Ulc ) and 
63 = El + 2 Ugj — 3J ex /2 + 2|£ c | + 16£^/ (Uq — Ulc ) un¬ 
der the condition Uo > Ulc where J ex = E 0 (S = 

3/2 )—Eo(S = 1/2) is the difference in energy between the 
spin-3/2 and spin-1/2 configuration. 13 However, the rel¬ 
ative strengths of peaks in the conductance spectrum can 
only be obtained by solving the full Anderson-Hubbard 
model self-consistently. 

Unlike the G e spectrum of DQDs,E 7 -l the G e spectrum 
of TQDM does not show the electron-hole symmetry due 
to topological effect. Note that N = 4 and N = 5 cor¬ 
respond to two-hole and one-hole configurations, respec¬ 
tively. A large Coulomb blockade separation between e3 
and 64 is given by A 34 = Uo + 3 J ex — 4t c — 8£^/(Uo — 
Ulc)- Here, 64 corresponds to the two hole ground state 
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with spin triplet instead of singlet. The magnitude of 
G e is smaller than the quantum conductance 2 e 2 /h for 
t LR /V > 1 as a result of electron Coulomb interactions J 17 
The mechanism for understanding the unusual G e be¬ 
havior in nanostructure junction systems is a subject of 
high interest W 1 Due to electron Coulomb interactions, 
the magnitudes of peaks are related to the probability 
weights of quantum paths, which are related to single¬ 
particle occupation numbers and many-particle correla¬ 
tion functions . 18 

To reveal the configurations for each main peak, the 
one-particle occupation number = (n^), inter¬ 
dot two particle correlation functions , and 

three particle correlation functions ((riL,a^c,-a^R,a}^ 
and (nL,a^c,-a^R,a}) are plotted in Fig. 3(b). We al¬ 
ways have the relation N R ^ = Nr i(T ^ because 

the outer QDs are directly coupled to electrodes, but not 
the central QD. Such a relation also holds for two-particle 
and three-particle correlation functions. The six main 
peaks in Fig. 3(a) indicate the filling of TQDM up to 
the n-electron ground state for n = I,-- - , 6 . For ex¬ 
ample, 62 indicates the formation of two-electron state 
with spin singlet, while 63 indicates the formation of 
three-particle state with total spin S = 1/2, which can 
be described as the spin-frustration stat e ! 1 1 l 1 2 i28J > Be- 
cause the on-site Coulomb interaction favors homoge¬ 
neous distribution of three electrons in TQDM, whereas 
the interdot Coulomb repulsion favors the charge fluc¬ 
tuation. As seen in Fig. 3(b), for < 63 ( 78r 0 ) JVg 
in each dot clearly displays the charge fluctuation be¬ 
havior. When TQDM goes into a three-particle state 
(A, > 63 ), the charge fluctuation is suppressed, and 
each QD is filled with one particle (with Nl j(T = N R ^ a = 
N c , a = 0.5), while (n Lia n c ,an R - a }^(n L ,anc,-an Ri<T ) = 
(nL,- a nc,a n R,-a }• This also demonstrates the spin frus¬ 
tration condition as depicted in the inset of Fig. 3(a). 

Figure 4 shows the charge stability diagram for zero- 
bias electrical conductance (G e ) and total occupation 
number (TV) as functions of gate voltages exerted on any 
two QDs (labeled by V g \ and V g 2 ) for a TQDM connected 
to three terminals. The magnitudes of G e and N are indi¬ 
cated by different colors. It is noticed that G e is enhanced 
on the borders that separate domains of different values 
of occupation number (TV) with larger G e occurring at 
V g i = V^2 • This is a result of higher degeneracy and 
charge-fluctuation in the state. The largest G e for N < 3 
occurs at the junction between N = 1 and N = 2 domains 
when V g i = V g 2 . This feature corresponds to the 62 peak 
of Fig. 3(a). The diagram Fig. 4(a) is simply a collection 
of curves displayed in Fig. 3(a) at different values of V g 2 
that shifts the QD energy levels. We note that in the 
domains of N = 1 and N = 2, the areas with stripes 


G e (e 2 /h) N 



v 9 i/r 0 v g 1 /r 0 


Fig. 4 Charge stability diagram of TQDM. Ft = To, 
ksT — 1.5r 0 , Ulc — Ucr — Ulr — 12Fo, Ej — 3Fo- The 
energy levels are shifted according to 

Ei — Ef — Ulr - X)m=i Pt,meV g m , where the gate coupling 
constants are /3l, 1 — 0.5, Pc, 1 = 1 , Pr, 1 = 0.5, /3l, 2 = 1 , 

Pc ,2 = 0.5, and p Ri 2 = 1. 


are not symmetrical with respect to gate voltage. In Ref. 
HD a capacitive interaction model was employed to plot 
the diagram of N. In their model, the electron hopping 
strength t^j was ignored. Consequently, the charge sta¬ 
bility diagram of G e cannot be obtained. The charge 
stability diagram of G e obtained by our full calculation 
[as shown in Fig. 4(a)] bears close resemblance to the 
experimental results as shown in Fig. 2 of Ref. Hi] 

So far, the results shown in Figs. 2-4 are all related 
to the linear response. To further clarify the QI effect 
at finite bias (V a = lOTo) for different temperatures, we 
plot in Fig. 5(a) the tunneling current as a function of 
center QD energy, Ec = Ep + A c for the configura¬ 
tion shown in the inset of Fig. 5(b). For A c > 2To, 
the tunneling current is suppressed as temperature in¬ 
creases. This is attributed to a reduction of electron pop¬ 
ulation in the electrodes for electrons with energy near 
Ep + lOTo- For Ac < 2To, we notice that 7V c?cr quickly 
jumps to 0.4, indicating that the central QD is filled 
with charge, which causes an interdot Coulomb block¬ 
ade for electrons entering the left QD. (See the reduc¬ 
tion of Nl,(j in Fig. 5(b)). This explains the sharp dip 
of G e for A c < 2To in Fig. 5(a). As temperature in¬ 
creases, such a dip in tunneling current is smeared out. 
At Ac = lOTo + t LC tcR/t L R = A qi, the tunneling cur¬ 
rent vanishes for all temperatures considered due to the 
QI effect. Such a robust destructive QI effect with respect 
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Fig. 5 (a) Tunneling current as a function of 
Ac = Ec — Ep at finite bias V a = 10To for various 
temperatures, (b) Occupation numbers as a function of Ac 
at ksT — 0. The other physical parameters used are 
El — Er — Er + 10To, Ulc = Ucr — 30To, Ulr — 10To, 
tLC — tcR — 3To, tLR — 0.4To, and T — To. 

to temperature provides a remarkable advantage for the 
realization of single electron QI transistors at room tem¬ 
perature. 12 To understand the interdot correlation effect, 
we also plot the case without interdot Coulomb interac¬ 
tion (Uij = 0) for ksT = 0 in Fig. 5(a). (See solid 
curve) We notice that the QI effect remains qualitatively 
the same, except that the tunneling current is slightly 
enhanced with interdot Coulomb interaction turned off. 
As the QI effect suppresses the current flow, the charge 
will accumulate in the left dot. Thus, N L ^ a reaches the 
maximum at Ac = Aqj while Nr reaches the minimum 
as seen in Fig. 5(b). This implies that the QI effect can 
be utilized to control charge storage in TQDM. 

4 Conclusions 

In summary we have obtained full solution to the charge 
transport through TQDM junction in the presence of 
electron Coulomb interactions, which includes all n- 
electron (n = 1 , • • • , 6 ) Green’s functions and correlation 
functions. The destructive and constructive QI behaviors 
of TQDM are clarified by considering the LDCT effect on 
the conductance spectrum. The conductance spectrum of 
TQDM with total occupation number varying from one 
to six directly reveals the electron-hole asymmetry due 
to topological effect. The calculated correlation func¬ 
tions also illustrate the charge fluctuation and spin frus¬ 
tration behaviors of TQDM. Our numerical results for 
charge stability diagram match experimental measure¬ 


ments very well. Finally, we demonstrated that the QI 
effect in TQDM is robust against temperature variation 
and it can be utilized to control the charge storage. 
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